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1 Introduction 



A major phenomenon that the theory of strong interactions, quantum chromodynamics 
(QCD), should account for is confinement: quarks and gluons are not observable parti- 
cles. In fact, every physical process involving strong interactions at high energy results 
in the formation of hadrons, in which quarks and gluons are confined on a distance scale 
of 0(1) fm. While, up to now, there is no formal proof that QCD implies confinement, 
there are many indications, both from perturbative and from lattice numerical studies, 
that this is likely the case. Perturbative QCD is applicable to scattering processes of 
quarks and gluons involving large momentum transfer (^> 1 GeV) because the strong 
coupling constant as is small enough to allow a series expansion. However, this is no 
longer possible at a scale of 1 GeV or below, where the perturbative expansion is mean- 
ingless, and where confinement and hadronization, the process of hadron formation, 
takes place. Thus, hadronization is not yet calculable from QCD first principles and 
one has to resort to phenomenological models. While this may seem an inconvenient 
limitation, still much can be learned from these models about QCD in the confinement 
regime. Indeed, if they are able to effectively describe the essential features of the ac- 
tual physical process, they give us relevant information about the characteristics of the 
fundamental theory. 

In this chapter, we will review two of these models, that have found widespread use in 
relativistic heavy ion collisions. The first is a model with a rather long history that has 
recently been revived by its successes in the description of hadronic multiplicities, the 
statistical model. This model is applicable to hadronization in elementary collisions as 
well as heavy ion collisions. This is in fact its main strength, i.e. it captures a universal 
feature in the hadronization process. 

The second model is the quark recombination or coalescence model which extends 
the concept of single parton fragmentation function, which has been used in elementary 
collisions since the '70s. Its recent success comes from observations specific to relativistic 
heavy ion collisions. 

We will start by reviewing the foundations and the main results of the statistical 
model in Sect. 2 and of the quark recombination model in Sect. 3. In Sect. 4 we will 
compare the two models and discuss further perspectives in the understanding of the 
hadronization phenomenon. 

2 The Statistical Hadronization Model 

The idea of applying statistical concepts to the problem of multi-particle production 
in high energy collisions dates back to a work of Fermi [1] in 1950, who assumed that 
particles originated from an excited region evenly occupying all available phase space 
states. This was one of Fermi's favorite ideas and soon led to an intense effort in trying 
to work out the predictions of inclusive particle rates calculating, analytically and nu- 
merically, the involved multidimensional phase-space integrals. When it became clear 
that the (quasi) isotropic particle emission in the center-of-mass frame predicted by 
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Fermi's model was ruled out by the data, an amendment was put forward by Hage- 
dorn [2] in the '60s, who postulated the existence of two hadron emitting sources flying 
apart longitudinally in the center-of-mass frame of a pp collision. Thereby, one could 
explain the striking difference between spectra in transverse and longitudinal momen- 
tum. Hagedorn was also able to explain the almost universal slope of pt spectra in his 
renowned statistical bootstrap model, assuming that resonances are made of hadrons 
and resonances in turn. 

After QCD turned up, many phenomenological models of strong interactions were no 
longer pursued and the statistical model was no exception. The resurgence of interest 
in these ideas came about when it was argued that a completely equilibrated hadron 
gas would be a clear signature of the formation of a transient Quark-Gluon Plasma 
(QGP) in heavy ion collisions at high energy. While it has been indeed confirmed that 
an (almost) fully equilibrated hadron gas has been produced [3] in those collisions, the 
interest in this model was also revived by the unexpected observation that it is able 
to accurately reproduce particle multiplicities in elementary collisions [I]. Naively, one 
did not expect a statistical approach to work in an environment where the number of 
particles is O(10) because it was a belief of many that a hadronic thermalization process 
would take a long time if driven by hadronic collisions. Apparently this is not the case 
and one of the burning questions, which is still waiting a generally accepted answer, is 
why a supposedly non-thermal system exhibits a striking thermal behavior. 



Figure 1: High energy collisions are assumed to give rise to multiple clusters at the 
hadronization stage [top]. Each cluster [bottom] is a colorless extended massive ob- 
ject endowed with abelian charges (electric, strange, baryonic etc.), intrinsic angular 
momentum and other quantum numbers such as parity, C-parity and isospin. 





Q = (Q,B,S....) 
M = mass 

J, X = spin, 3 ,d comp. 

V = volume 

Isospin, parity, C-paritv 
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Before we address this interesting issue, it is appropriate to provide a rigorous for- 
mulation of the model in a modern form, which is necessarily different from Fermi's 
original model due to the tremendous improvement in our knowledge of strong interac- 
tions phenomenology. The Statistical Hadronization Model (SHM) must be considered 
as an effective model describing the process of hadron formation in high energy collisions 
at energy (or distance) scales where perturbative QCD is no longer applicable. A high 
energy collision is thought of as a complex dynamical process, governed by QCD, which 
eventually gives rise to the formation of extended massive colorless objects defined as 
clusters or fireballs (see Fig. [I]). While the multiplicity, masses, momenta and charges 
of these objects are determined by this complex dynamical process, the SHM postulates 
that hadrons are formed from the decay of each cluster in a purely statistical fashion, 
that is: 

Every multihadronic state localized within the cluster and compatible with 
conservation laws is equally likely. 

This is the urprinzip of the SHM. The assumption of the eventual formation of 
massive colorless clusters is common for many hadronization models (e.g. the cluster 
model implemented in the Monte-Carlo code HERWIG [5]) based on the property of 
color preconfinement [Hj exhibited by perturbative QCD. The distinctive feature of the 
SHM is that clusters have a finite spacial size. This aspect of clusters as a relativistic 
massive extended objects coincides with that of a bag in the MIT bag model [7]. Indeed, 
the SHM can be considered as an effective model to calculate bag decays. 

The requirement of finite spacial extension is crucial. If the SHM is to be an effective 
representation of the QCD-driven dynamical hadronization process, this characteristic 
must be ultimately related to the QCD fundamental scale Aqcd- As we will see, the 
universal soft scale shows up in the approximately constant energy density at hadroniza- 
tion; in other words, the volume of clusters is in a constant ratio with their mass when 
hadronization takes place. It is also worth stressing here that there is clear, indepen- 
dent evidence of the finite size of hadronic sources in high energy collisions. Quantum 
interference effects in the production of identical particles, the so-called Bose-Einstein 
correlations or Hanbury Brown-Twiss second-order interference, is by now a firmly es- 
tablished phenomenon. This effect would simply be impossible without a finite volume. 

2.1 Localized States 

The basic postulate of the Statistical Hadronization Model asserts that every localized 
multihadronic state which is contained within a cluster and is compatible with conser- 
vation laws is equally likely. The word localized, implying a finite spacial size, plays 
a crucial role, as we have emphasized. Thus, before getting to the heart of the SHM 
formalism, it is necessary to pause and clarify the distinction between localized and 
asymptotic states. 

Such a difference is not an issue when the volume is sufficiently larger than the 
Compton wavelength of hadrons and it is disregarded in most applications where clusters 
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supposedly meet this requirement (e.g. heavy ion collisions); yet, it is an important 
point at a fundamental level. Although in thermodynamics the focus is on the limit of 
infinite volumes, we must start from a finite volume and localized systems to introduce 
concepts like energy density, temperature etc. Furthermore, in the hadronic world, finite 
size effects must be diligently taken into account when the volume is comparable to the 
(third power of) the pion's Compton wavelength, ~ 1.4 fm. 



Figure 2: The localized multi-hadronic states \hy) pertaining to the quantum field 
problem in a limited region. Asymptotic states |/) are the usual free states characterized 
by particle momenta and spin components. 

The difference between a localized and an asymptotic state is depicted in Fig. [2j 
For a single particle in a Non-Relativistic Quantum Mechanics (NRQM) framework, the 
conceptual difference is easier to grasp: a localized state is described by a wavefunction 
which vanishes outside the cluster's region whereas an asymptotic state is a wavefunction 
which is defined over the whole space (e.g. a plane or a spherical wave). In Quantum 
Field Theory (QFT), a localized state is a state of the Hilbert space defined by the 
localized problem, e.g. the problem of the quantum field in the cluster's finite region. 
For a free field, if we enforce fixed or periodic boundary conditions, such states are simply 
defined by integer occupation numbers for each allowed mode in the finite region, as is 
well known. For a multiparticle state of non-interacting hadrons, this state will be 
defined by all occupation numbers of the modes determined by the fields associated 
with the different species of hadrons, and we will simply denote it with \hy) (where h 
stands for "hadrons" and V stands for the finite region of volume V). On the other 
hand, an asymptotic state is a state of the Hilbert space defined by the quantum field 
operators over the whole space; for a free field, these are the familiar multi-particle free 
states defined by, e.g., momentum and polarization: 



There is a noteworthy and deep difference between non- relativist ic quantum mechanical 
and quantum field theoretical case. In the latter, particle number is not fixed and 
a localized state with multiplicity (defined as the sum of all occupation numbers) N 




\f) = \Pi,<?i, ■ ■ .,Pn,(?n) 
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does not necessarily correspond to an asymptotic state with particle multiplicity A^|. 
Unlike in NRQM, a localized state with multiplicity N has non-vanishing projections 
over asymptotic states with different particle multiplicities. In symbols: 

\N) V = a 0tN \0) + a hN \l) + ... + a N , N \N) + ... (2) 

with the obvious condition that a^jy — ► when V — > oo for i ^ N . In particular, the 
vacuum of the finite-region problem |0)y is different from the vacuum of the full-space 
problem |0), which is commonly known as the Casimir effect. With a straightforward 
mapping of the Hilbert space of the localized quantum field onto the full Hilbert space, 
it is possible to express destruction operators of the localized field as linear combinations 
of destruction and creation operators of the field defined over the whole space [8]. These 
relations would be sufficient to calculate the coefficients of the above equation, but this 
is not really needed in the SHM, as it will be soon clear. 



2.2 The Formalism: Basics 

Let us consider a cluster and assume first that it can be described as a mixture of states. 
Then, the basic postulate implies that the corresponding density matrix is a sum over 
all localized states projected onto the initial cluster's quantum numbers: 

poc ^Pi\h v ){h v \Pi = PiPyPi (3) 

hv 

where \hy) are multi-hadronic localized states and P« is the projector onto the cluster's 
initial conserved quantities: energy-momentum, intrinsic angular momentum and its 
third component, parity and the generators of inner symmetries of strong interactions 

a 

The operator Pj can be formally defined as the projector onto an irreducible vector of 
the full symmetry group and worked out in a group theory framework [HJ [TU| . It can be 
factorized into a "kinematic" projector, associated to general space-time symmetry, and 
a projector for inner symmetries. For the space-time symmetry, the relevant group is the 
extended orthochronous Poincare group 10(1,3)^ and an irreducible state is defined by 
a four-momentum P, a spin J and its third component A and a discrete parity quantum 
number 7T = ±1. Therefore: 

Pj P P,J,X,nP inner (4) 

If the projector Pp,j,\ is worked out in the cluster's rest frame where P = (M, 0), it 
further factorizes into the product of simpler projectors 0,(11], i.e.: 

Pp,J,X,, = S\P-P)Pj, X ^^ (5) 

1 We note that the multiplicity of an asymptotic free state is the properly defined number of particles 
2 Opcrators in the Hilbert space will be denoted with a hat. Exceptions to this rule are projectors, 
which will be written in serif font. i.e. P. 
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where P is the four-momentum operator, Pj\ is a projector onto SU(2) irreducible states 
| J, A) and II is the space reflection operator. 

As clusters are color singlets by definition, the projector Pi nner involves flavor and 
baryon number conservation. In principle, the largest symmetry group one should con- 
sider is SU(3), plus three other U(l) groups for baryon number, charm and beauty 
conservation. However, SU(3) symmetry is badly broken by the mass difference be- 
tween strange and up, down quarks, so it is customary to take a reduced SU(2)®U(1) 
where SU(2) is associated with isospin and U(l) with strangeness. The isospin SU(2) 
symmetry is explicitly broken as well, but the breaking term is small and can generally 
be neglected. However, most calculations in the past have replaced isospin SU(2) with 
another U(l) group for electric charge, so that the symmetry scheme, from an original 

SU(2)j SOS pj n ®U(l) s £rangeness®U(l)h ar y 0n reduces to U (1 .) charge®^ (-0 strangeness®^ ^ ( 1 ) baryon • 

Altogether, Pi nner can be written as 

Pinner = P/,/ 3 PqPx (6) 

where / and Is are isospin and its third component, Q = (Qi, ■ ■ ■ ,Qm) is a vector of 
M integer abelian charges (baryon number, strangeness, etc.) and P x is the projector 
onto C-parity, which makes sense only if the system is completely neutral, i.e. 1 = 
and Q = 0; in this case, P x commutes with all other projectors. 

From the density matrix ([3]) the probability of observing an asymptotic multiparticle 
state \f) is 

P/<x</|P,PyP,|/> (7) 

which is well-defined in terms of positivity and conservation laws. In fact, ((7j) is mani- 
festly positive definite and pf = if the state |/) has not the same quantum numbers as 
the initial state. By summing over all states |/), one obtains the trace of the operator 
PiPyPj which is 

J^Pf « tr^PyP;) = tr(P 4 2 P y ) = atr(P i P y ) . (8) 

/ 

The constant a is divergent and positive. It can be directly checked by choosing the |/) 
as momentum eigenstates and using the expression on the right hand side of (j3J). The 
reason for its presence is the non-compactness of the Poincare group, which makes it 
impossible to have a properly normalized projector. The last trace in (jSJ) can be written 

as 

tT(p i p v ) = ^2{h v \p i \h v ) = n (9) 

hv 

which is, by definition the microcanonical partition function [5], i.e. the sum over all 
localized states projected onto the conserved quantities defined by the selected initial 
state. If only energy and momentum conservation is enforced, Q takes on a more familiar 
form: 

tt = ^2(h v \5 4 (P-P)\h v ). (10) 

hv 
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Although the mixture of states defined by Eq. allows us to calculate probabilities 
of any measurement unambiguously, a cluster could in principle also be described by a 
pure quantum state. Actually, the mixture of states only expresses our ignorance about 
the true state of the system, which is in principle a pure one, or, more precisely, a pure 
state entangled with pure states of other clusters. To avoid slipping into fundamental 
quantum mechanics problems of decoherence and measurement, we take a pragmatic 
stance here. It suffices to realize that in some low-energy collision events, only one 
cluster might be created whose state is then necessarily a pure one. According to the 
postulate of the SHM, this must be an even superposition of all localized states with 
the initial conserved quantities, i.e. 

\^P) = S ^c hv Pi\h v ) with \c hv | 2 = const . (11) 

h \- 

The probability of observing a final state \f) is then 

iw>i 2 = iE</i p *i^> c *vi a (12) 

hy 

= const^|(/|P,|M| 2 + E (/I^IM^ylPd/)^^- 

hy hv^h'y 

If the coefficients Ch v have random phases, the last term in Eq. (|12p vanishes and we are 
left with the same expression appearing in Eq. (J7J); in other words an effective mixture 
description is recovered. Hence a new hypothesis is introduced in the SHM here: if the 
cluster is a pure state, the superposition of multi-hadronic localized states must have 
random phases. 

Now the main goal of the model is to determine the probabilities fl7]) which involves 
the calculation of the projector Py = J2h v \hv)(hv\, a more limited task than the 
explicit calculation of all scalar products (hy\f). Since the states \hy) are a complete 
set of states of the Hilbert space Hy for the localized problem, the above projector is 
simply a resolution of the identity of the localized problem and can be written in the 
basis of the field states. For a real scalar field this is 

Py= f £>V#)(V>I (13) 

Jv 

where = ® X |V'( X )) and T>ip is the functional measure; the index V means that the 
functional integration must be performed over the field degrees of freedom in the region 
V, that is Dip = rixGy^( x )- One nas ^° l ace severa l conceptual subtleties in the 
endeavor of calculating the probabilities ((7j) with the projector (fT3|) . e.g. how to deal 
with field boundary conditions and with their values outside the region V. However, by 
enforcing the known non-relativistic limit is possible to come to an unambiguous and 
consistent result [8]. 
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2.3 Rates of Multiparticle Channels 



According to the formulae introduced in the previous section, the decay rate of a cluster 
into a channel {Nj} ({Nj} is the array of multiplicities (Ni, . . . , Nk) for hadron species 
1,...,K) is proportional to the right hand side of (jTj) integrated over momenta and 
summed over polarization states of the final hadrons. Taking into account only energy- 
momentum conservation, so that the projector (jSJ) reduces to 

Ppjx, - Pp = 5 4 (P - P) , 

and neglecting quantum statistics effects, this is proportional to the microcanonical 
partition function with fixed particle multiplicities [HI [8] 



n 




d 3 Pl ... j d 3 p N 5\P -J2p*)(0\Pv\0). 



(14) 



Here N is the number of particles, Sj the spin, P = (M, 0), M is the mass and V is the 
cluster's proper volume. This formula is the same as it would be obtained in NRQM, 
with the factor (0|Py|0) (which becomes 1 in the limit V — > oo) being the only effect 
of the field theoretical treatment [8]. Since only relative rates make sense, this common 
factor for all channels is irrelevant. 

Loosely speaking, Eq. (1141) tells us that the decay rate of a massive cluster into some 
multi-hadronic channel is proportional to its phase space volume. However, it should 
be emphasized that the "phase space volume" in ()14p is calculated with the measure 
d 3 xd 3 p for each particle, and not with the one usually understood in QFT, i.e. d 3 p/2e. 
Although this is also commonly known as "phase space", it is quantitatively different 
from the properly called phase space measure d 3 x d 3 p and should be called " invariant 
momentum space" measure |12j . 

Eq. (fT4"]) can be cast in a form which makes its Lorentz invariance apparent. Define 
a four-volume T = Vu [12] where V is the cluster's rest frame and u its four-velocity 
vector. Then (1141) can be rewritten as: 



{Nj} 



1 



(2yr) 3iV 

N 



n 



(2Sj + 1 
NA 



\T ■ Vi 8{pt - m1)e(j>\ 




i=l 



d A pi ... j d 4 p N 



5 4 (P -$>)(0|Py|0> 



(15) 



which is manifestly covariant. In this form it can be directly compared with the general 
formula for the decay rate of a massive particle into a iV-body channel: 



.,(7 AT 




d 3 p N 
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N 



\M„\H S (P„ 



E 



Pi 



(16) 
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where a labels, as usual, polarization states. Comparing (|14p with (fTol) we can infer a 
dynamical matrix element for the SHM which is 



N 



|M /4 | 2 oc l[T Pl (17) 



i=l 



Therefore, according to the SHM the dynamics in cluster decay is limited to a common 
factor for each emitted particle, which linearly depends on the cluster's spacial size. The 
four-volume T is simply proportional to the four-momentum of the cluster through the 
inverse of energy density p and therefore: 



1 N 

JtU P ^ ( 18 ) 



I f%\ p 



This expression explicitly shows the separation between the kinematic arguments of the 
dynamical matrix element, and the scale 1/p which determines particle production. This 
ought to be ultimately related to the fundamental scale of quantum chromodynamics, 
Aqcd- 

It has already been stressed that the finite cluster size is the distinctive feature of 
the SHM. This peculiarity of the model stands out when taking into account quantum 
statistics for the calculation of decay rates. Our final result, for which (1141) is a special 
case when all particles belong to different species, reads 

= / dv..rfwvo-i>)nE (T1 'T ( s, + '' Hj (19 » 

J i=l 3 {h nj } [[n 3 =l n j S ! 

Hi 

x 11^ <0|Rv-|0> 

lj=l 

where {h n .} is a partition of the integer Nj in the multiplicity representation, that is 

E5 i " A , = iN 

Fourier integrals: 

F nt = J] d 3 x jxiP^o-Pn) (20) 

over the cluster's region V, q being a cyclic permutation of order n\. The expression 
( TL9]) has been obtained in refs. [91 [8] and is a generalization of a similar one calculated 
by Chaichian, Hagedorn and Hayashi [T2] whose validity is restricted to large volumes. 
It is a so-called cluster decomposition of the microcanonical partition function of the 
channel. For sufficiently large volumes, all terms in Eq. f|l~9l) turn out to be proportional 
to the Hjth power of the volume V [9] , so that the leading term is the one with Hj = Nj 
for all j, which leads precisely to Eq. (JHj) . Thus, Eq. f|l~9l) is a generalization of ([Tj 
containing all corrections due to quantum statistics. 



X];); i ",/',., = Nj, J2n-=ihnj = Hj and Y.j N i = N - Tne factor F ni , in Eq. §F5§ are 
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In general, with respect to the Boltzmann case (|14p . the channel rate is enhanced 
by the presence of identical bosons and suppressed by that of fermions. This means 
that Bose-Einstein and Fermi-Dirac correlations are built into the SHM. The reader has 
probably anticipated this fact through the appearance of typical Fourier integrals in the 
cluster decomposition. This feature of the SHM is an almost obvious consequence of the 
cluster's finite spacial size. 

2.4 Interactions 

So far, we have dealt with non-interacting particles. However, the localized hadronic 
fields we have used to calculate transition probabilities do interact and this must be 
taken into account. The energy of the interacting system must be conserved until the 
final asymptotic multi-hadronic state is reached which is made of particles stable under 
strong interactions, namely pions, kaons, nucleons and octet hyperons. 

Formally, this implies that the projector (jSJ) must include the interacting Hamil- 
tonian in the <5 4 (P — P) operatoiEl The definition (JTj) for the probability to observe 
an asymptotic state |/), has to be modified by the insertion of M0ller operator Q, yet 
summing over the complete set of states yields the same result as in Eq. (jSJ): 



where a is an irrelevant divergent constant and Q is the microcanonical partition function 
of the interacting hadronic system. 

An outstanding theorem by Dashen, Ma and Bernstein (DMB) [T3] allows us to 
calculate the microcanonical partition function of an interacting system in the thermo- 
dynamic limit V — > oo as the sum of the free one plus a term depending only on the 
physical scattering matrix. It can be expressed as 



where P includes the full interaction Hamiltonian, whereas Pq only contains the free one; 
S is the reduced scattering matrix on the energy-momentum shell. If more conserved 
quantities other than energy and momentum are involved, like those encountered in 
section [2~2l the theorem is readily extended and relevant projectors can be placed next 
to the (^-functions in Eq. (1221 ; it suffices that these conserved quantities are associated 
with symmetries of both free and interacting theory. 

This theorem is indeed the starting point of the hadron-resonance gas model since 
it can be shown that if only the resonant part of the scattering matrix is retained and 
the background interaction can be neglected, the main contribution of the second term 

3 In all virtually known field theories, there is no additional interacting term for momentum and 
angular momentum |13| . 




(21) 




(22) 
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Figure 3: Left panel: symmetric diagrams for the cluster decomposition of the interaction 
term in the DMB theorem. Right panel: non-symmetric diagrams. 

on the right hand side of Eq. (1221) is equivalent to considering all hadronic resonances 
as free particles with distributed mass. More specifically, if a cluster decomposition of 
the two scattering operators is carried out in Eq. ( J22l) . the corresponding diagrams can 
be divided in two sets: the symmetric diagrams (see Fig. [31 left panel) and the non- 
symmetric ones (see Fig. [3l right panel). Taking into account that the terminal legs on 
both sides have to be the same stable particles on entry and exit (we are calculating 
a trace), it can be shown that the main contribution to symmetric diagrams comes 
from the matching resonances in bubbles facing each other. For each term of the trace, 
this amounts to adding the decay products of resonances considered as free particles 
with masses distributed according to a relativistic Breit-Wigner form. In symmetric 
diagrams, there is in principle an additional contribution from resonance interference, 
which might be non-negligible in case of wide, overlapping resonances with the same 
decay channel, but it depends on mostly unknown complex parameters and it is thus 
neglected. 

Likewise, the asymmetric diagrams give an additional contribution which also de- 
pends on the aforementioned complex interference parameters. While the number of 
such diagrams greatly exceeds the symmetric ones due to the large number of reso- 
nances, contributing terms can be both positive and negative and hopefully a partial 
cancellation occurs when summing them up for a selected final state. 

Altogether, retaining only the resonant interaction and symmetric diagrams in the 
scattering matrix cluster decomposition, and neglecting resonance interference leads to 
the following picture: an interacting hadron gas is, to a good approximation, a gas of 
non-interacting free hadrons and resonances. Since non-resonant interaction should be 
negligible, the ideal hadron-resonance gas picture holds if the energy density or temper- 
ature of the system is large enough for most resonances to be excited. A quantitative 
assessment of how large these parameters are is still missing, a rough estimate being 
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T > 100 MeV. 

An important remark is now in order. The DMB theorem affirms the equality of two 
traces, but not of single trace terms. Yet, the decomposition of Eq. (!22j) . implying the 
ideal hadron-resonance gas picture, is widely used for the calculation of inclusive stable 
hadronic multiplicities as well, which requires a condition stronger than the equality 
of the traces on both sides. In other words, using the decomposition ( 1221) to calculate 
average multiplicities or fluctuations requires the equality to hold for multiparticle gen- 
erating functions and not only for microcanonical partition functions. Up to now, the 
extension of fl22l) to generating functions has never been proved; most likely, it does not 
hold and corrections to this assumption are necessary. Moreover, while the theorem re- 
quires the thermodynamic limit, it is commonly used at finite volume. These limitations 
should be always kept in mind when using the ideal hadron-resonance gas model. 

2.5 High Energy Collisions 

As we have seen in Sect. 12.31 each individual cluster produced in a high energy collision 
(shown in Fig. [1]), should be hadronized according to formula ( TT9l) . or its approximation 
(TH|) . which yields the rates of a given channel within the microcanonical ensemble, in- 
cluding energy-momentum conservation. If clusters are large enough, the microcanonical 
ensemble could be well approximated by the canonical [HI [15] or even grand-canonical 
ensemble for average multiplicities. This is not the case in elementary collisions (e + e~, 
pp, etc.) while it is generally possible in heavy ion collisions, as we will see later in this 
section. Calculating observables in high energy collisions within the SHM then implies 
summing microcanonical averages over all produced clusters and this requires in turn 
knowledge of their charges and four-momenta. In fact, this latter information is un- 
known to the SHM and only a dynamical model of the pre-hadronization stage of the 
process (such as, e.g., HERWIG) can provide it. 

However, if we are interested in calculating Lorentz-invariant observables (such as 
average multiplicities) the momenta of clusters are immaterial and only charges and 
masses matter. In this case, one can introduce a peculiar extra-assumption which allows 
to considerably simplify the calculation. Basically, it is assumed that the probability 
distribution: 

w(Q 1 ,M 1 ;...,Q N ,M N ) 

of masses M and conserved abelian charges Q for N different clusters is the same 
as one would have by randomly splitting a large cluster (defined as Equivalent Global 
Cluster, EGC) into N subsystems with given volumes. Thereby, the Lorentz invariant 
observables can be calculated for one (equivalent global) cluster, whose volume is the 
sum of proper cluster volumes and whose charge is the sum of cluster charges, hence 
the conserved charge of the initial colliding system. The full mathematical procedure is 
described in detail in ref. [T6] . 

In such a global averaging process, the EGC generally turns out to be large enough 
in mass and volume so that the canonical ensemble becomes a good approximation 
of the more fundamental microcanonical ensemble [15]; in other words, a temperature 
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can be introduced which replaces the a priori more fundamental description in terms 
of energy density. It was shown that the mass of the EGC should be at least 8 GeV 
(with an energy density of 0.5 GeV/fm 3 ) for the canonical ensemble to be a reasonably 
good approximation [15] . Also, it should be emphasized that in such a mathematical 
reduction process, temperature has essentially a global meaning and not local as in 
hydrodynamical models (see next subsection). The only meaningful local quantity in 
actual physical process are energy densities and individual physical clusters cannot be 
described in terms of a temperature, unless they are sufficiently large. Nevertheless, 
this "global" temperature closely mirrors the value of energy density at which clusters 
hadronize. Indeed, it is this latter value which mainly determines hadronization-related 
observables; the requirement of the charge distribution of EGC is a side-assumption 
which is important to simplify calculations, but it can possibly be replaced by other 
distributions leaving final results essentially unchanged. 

In this approach, the primary multiplicity of each hadron species j is given by [16]: 

/„ \P~y _ VT(2Sj + 1) ^ „ +1 m) / n mA ZjQ-nqj) 



n=l 



where V is the (mean) volume and T the temperature of the equivalent global clus- 
ter. Here Z(Q) is the canonical partition function depending on the initial abelian 
charges Q = (Q, N, S,C, B), i.e., electric charge, baryon number, strangeness, charm 
and beauty, respectively; rrij and Sj are the mass and the spin of the hadron j, q^- = 
(Qj, Nj, Sj, Cj, Bj) its corresponding charges; the upper sign applies to bosons and the 
lower sign to fermions. 

The parameter 75 in (123]) is an extra phenomenological factor implementing an ad hoc 
suppression of hadrons with N s strange valence quarks with respect to the equilibrium 
value. This parameter is outside a pure thermodynamical framework and it is needed 
to reproduce the data, as we will see. For temperature values of 160 MeV or higher, 
Boltzmann statistics, corresponding to the term n — 1 only in the series (l2"3l . is a very 
good approximation (within 1.5%) for all hadrons but pions. For resonances, the formula 
( 123]) is folded with a relativistic Breit-Wigner distribution for the mass rrij. The final 
multiplicities, to be compared with the data, are determined by adding to the primary 
multiplicity (123]) the contribution from the decay of unstable heavier hadrons, according 
to the formula 

(n,) = (n t y~y + J2Mk^j)(n k ) . (24) 

k 

The canonical partition function can be expressed as a multi-dimensional integral 



x exp 



V 
(2 



V — Y^Sj + 1) J d 3 p log (1 ± ^e-V^m-i^+i 



(25) 
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Figure 4: Upper panel: measured vs theoretical multiplicities of light-flavoured hadrons 
in e + e~ collisions at y/s = 91.25 GeV. Lower panel: fit residuals (from ref. [T7]). 



where N is the number of conserved abelian charges. Unlike the grand-canonical case, 
the logarithm of the canonical partition function does not scale linearly with the volume. 
Therefore, the so-called chemical factors Z (Q — nqj) / Z (Q) [18] turn out to be less than 
unity even for a completely neutral system at finite volume (canonical suppression) and 
reach their grand-canonical value 1 at asymptotically large volumes [191 H] • 

The light-flavoured multiplicities in e + e~show a very good agreement with the pre- 
dictions of the model, as it shown in Fig. HI the temperature value is about 160 MeV 
and the strangeness undersaturation parameter 75 ~ 0.7. Similar good agreements are 
found for many kinds of high energy elementary collisions over a large energy range [1]. 
Also, an excellent agreement between measured and predicted relative abundances of 
heavy flavoured hadronic species in e + e~collisions by using the model parameters fitted 
to light-flavoured multiplicities (HHTj, as shown in Tabled! 

The overall striking feature is that the temperature turns out to be approximately 
constant over two orders of magnitude in centre-of-mass energy with a value of 160- 
170 MeV (see Fig. [5]) and very close to the QCD critical temperature as determined 
from lattice calculations. There must certainly be a profound connection between the 
thus-found hadronization temperature and QCD thermodynamics, a connection which 
has not been made clear yet. Nevertheless, this finding indicates that hadronization is 
a universal process occurring at a critical value of the local energy density, i.e. when 
clusters have an energy density of ~ 0.5 GeV/fm 3 . 

The parameter 75 is found to be less than 1 in all examined elementary collisions, 
ranging from ~ 0.5 in hadronic collisions to ~ 0.7 in e + e~collisions (see Fig. [9]). This 
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Figure 5: Temperatures fitted in elementary collisions as a function of center-of-mass 
energy. 



extra parameter most likely reflects the different mass of the strange quark with respect 
to lighter u, d quarks. This is a second scale, besides Aqcd, which must play a role 
in hadronization in view of its value (9(100) MeV. Altogether, one can say that the 
SHM description of hadronization is in excellent agreement with QCD at least with 
regard to the the number of parameters. The two parameters T and 75 correspond to 
the two fundamental scales Aqcd and m s , the strange quark mass. While we lack a 
definite relation connecting them (see, however ref. [17] ). it is worth stressing that a 
phenomenological description of hadronization in terms of fewer parameters cannot be 
possible. 

Finally, the statistical model shows a very good capability of reproducing transverse 
momentum spectra in hadronic [IE] as well as heavy ion collisions [20]. Particularly, the 
phenomenon of approximate mj> scaling observed in pp collisions is nicely accounted 
for by the model. However, the exact p-r conservation at low energy and the increasing 
importance of jet emission at high energy restrict the validity of the statistical canonical 
formulae to a limited centre-of-mass energy range. Within this region, a clear consistency 
is found between the temperature parameter extracted from the spectra and that from 
average multiplicities. Altogether, this finding bears out one of the key predictions of 
the SHM, namely the existence of a definite relation between the dependence of particle 
production rates on mass and, for each particle species, their momentum spectra (in 
the cluster's rest frame) because they are both governed by one parameter, the energy 
density (or temperature) at hadronization. 
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Table 1: Abundances of charmed hadrons in e + e~ — > cc annihilations and bottomed 
hadrons in e + e~ — > bb annihilations at y/s = 91.25 GeV, compared to the prediction of 
the statistical model (from ref. [T7]). 



2.6 Heavy Ion Collisions 

In heavy ion collisions, the system is much larger and two possibilities are usually en- 
visaged: either hadronizing clusters are simply much larger than those in elementary 
collisions; or clusters are hydrodynamical cells, i.e. they are small but in thermal con- 
tact with each other due to previous thermalization, which implies a strong correlation 
between their position and momentum and charge densities (see Fig. [6]). In both case 
the canonical or grand-canonical formalisms apply to individual clusters. For the former 
case, it is worth mentioning that the transition from a canonical to a grand-canonical 
description effectively occurs when the cluster volume is of the order of 100 fm 3 at an 
energy density of 0.5 GeV/fm 3 [TS]. If the EGC reduction assumption still applies, 
chemical factors in Eq. ([23"]) are replaced by fugacities, and in this case the phase-space 
integrated multiplicities read 

,2 



{n . r ^ = vt(2S, + d £ 7? ,„ (Tir+1 !2, K2 ^ exp[np . ^ (26) 

n=l 
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Figure 6: Spacial distributions of clusters in heavy ion collisions according to the hy- 
drodynamical picture. In this model, nearby clusters interact from an early stage on 
and their momenta and charges are strongly correlated with their positions, unlike in 
elementary collisions. 



fi is a vector of chemical potentials pertaining to the conserved abelian charges, i.e. the 
electrical chemical potential /xq, the baryon chemical potential /xg and the strangeness 
chemical potential /x^. Usually, but not always, /xg and /xq are determined by enforcing 
strangeness neutrality and by fixing the ratio Q/B to be the same as the initial Z/A 
ratio of the colliding nuclei. 

In the framework of the hydrodynamical model, formula (1261) applies to individual 
clusters identified with hydrodynamic cells and both temperature and chemical poten- 
tials depend on space-time; when integrating particle densities to get average multiplic- 
ities, one should take into account this dependence. It is important to stress that the 
hydrodynamical description is a salient feature of heavy ion collisions due to the early 
thermalization of the system in the partonic phase, a phenomenon which does not oc- 
cur in elementary collisions. It is this early thermalization which establishes the strong 
correlation between positions and velocities of clusters, supposedly absent in elementary 
collisions. 

Provided that rapidity distributions are wide enough, and that there is little variation 
of the thermodynamical parameters of clusters around midrapidity, the formula (T2B1) 
describes rapidity densities of hadrons at midrapidity as well: this condition is fulfilled 
at RHIC energies, but not at AGS and SPS energies, where the measured rapidity 
distributions are not significantly wider than those of a single fireball at the temperature 
found [22]. 

In general, the fits to particle multiplicities in heavy ion collisions are of the same 
good quality as in elementary collisions (see Fig. [7]). Many groups have analyzed the 
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Figure 7: Upper panel: measured vs theoretical multiplicities of light-flavoured hadrons 
in Pb-Pb collisions at y/s NN = 17.2 GeV. Lower panel: fit residuals (from ref. [2~T]). 



data over more than a decade [3] and the overall description is very good throughout 
all explored energies and one finds a smooth curve in the T — \i& plane (see Fig. [8]). 

2.7 Strangeness Production 

The statistical model is a very useful tool to study one of the main features of relativis- 
tic heavy-ion collisions, the increase of relative strangeness production with respect to 
elementary collisions. This was one of the early signatures proposed for Quark-Gluon 
Plasma formation [24], and it has therefor attracted much attention both on the theoret- 
ical and experimental side. According to the SHM, this is mainly an effect of the increase 
of the global volume from elementary to heavy-ion collisions. In elementary collisions, 
the EGG volume is small enough for the chemical factors (see Eq. ff25l) ) of strange 
particles to be consistently less than 1 for systems with vanishing net strangeness, a 
phenomenon known as strangeness canonical suppression. 

However, canonical suppression is not enough to account for strangeness enhance- 
ment from pp to heavy-ion collisions: also an increase of 75 is needed. This is demon- 
strated by neutral mesons containing strange quarks, especially <fi meson, which do not 
suffer canonical suppression but are relatively more abundant in heavy-ion collisions 
[25l [26] . Therefore, from a SHM viewpoint, one can say that, as far as particle abun- 
dances is concerned, the only substantial difference between elementary and heavy ion 
collisions resides in the different 75 values, which are generally higher in heavy-ion col- 
lisions and increase slowly as a function of center-of-mass energy (see Fig. [9]): at RHIC 
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energies one finds 75 ~ 1 in central collisions. However, since 75 is an empirical param- 
eter which lies outside of a pure statistical mechanics framework, this observation does 
not clarify the origin of strangeness enhancement. 

It is interesting to note that, while 7s shows no special regularity in elementary 
collisions, the ratio of newly produced ss pairs over one half uu + dd pairs (the so- 
called Wroblewski ratio \$) turns out to be around 0.2-0.25 at all energies, whereas it 
is definitely higher in heavy ion collisions (see Fig. ITUT) . 

These differences, both in 75- and Xs have been and still are subject of investigation. 
Relevant information comes from the centrality dependence of strangeness production, 
which provides an interpolation from single pp collisions at large values of nuclear impact 
parameter to head-on heavy-ion collisions at low values thereof. The enhancement has 
been measured by the experiments WA97 and NA57 at SPS energy [27J and STAR 
at RHIC [2B] for hyperons and other strange particles and it has been found to be 
hierarchical in strangeness content (highest for Q~ , lowest for A). These observations 
led some authors [29] to put forward a picture where 75 is an effective parametrization 
of a canonical suppression. For large enough baryon number and charge, it is possible to 
take a mixed canonical-grand-canonical approach where only strangeness conservation 
is enforced, while electric and baryon-chemical potential are introduced. The chemical 
factors Z(S—Sj)/Z(S) depend on the volume and saturate at large volumes, as expected. 
Therefore, if we want to account for 75 < 1 with this mechanism, there must be some 
small sub-regions within a large fireball where strangeness is exactly vanishing even for 
the most central collisions. Thereby, chemical factors are significantly less than 1 and a 
suppression with respect to the grand-canonical limit is implied. However, this model 
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Figure 9: Left panel: the strangeness undersaturation parameter 75 as a function of 
center-of-mass energy in e + e~, pp and pp collisions. Right panel: the strangeness un- 
dersaturation parameter 75 as a function of center-of-mass energy in central heavy ion 
collisions (from ref. |23j). 



has two major problems: 

1 . Since measured enhancement steadily increases from peripheral to central collisions 
and hadronization temperature does not change [301 EH E3], the volume of the 
sub-regions with S = should also increase and a saturation is thus expected (see 
Fig. [TTj) : yet, no saturation is observed, which is quite an oddity. 

2. As has been mentioned, canonical suppression has no effect on <p\ yet, the relative 
yield of this meson with two constituent strange quarks is also observed to increase 
from peripheral to central collisions [32J and with 75 = 1 and the observed constant 
temperature, this cannot occur. 

Recently, a geometrical explanation of these two features has been advocated [33] 
based on a superposition of emission from a hadron-resonance gas at full chemical equi- 
librium with 75 = 1, defined as the core, and from nucleon-nucleon collisions at the 
boundary of the overlapping region of the two colliding nuclei, defined as the corona, 
from which produced particles escape unscathed. Since in NN collisions strangeness is 
suppressed with respect to a fully equilibrated, grand-canonical hadron gas, if such NN 
collisions account for a significant fraction of total particle production, a global fit to 
particle multiplicities will find 75- < 1, as indeed observed in data. The idea of superpos- 
ing different sources is common to other models (a similar one is discussed in ref. |5^]). 
The peculiar feature of this specific model is to assume single NN collisions as secondary 
sources; only in this case does it seem possible to reproduce centrality dependence of 
the (j) meson. 
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Figure 10: Wroblewski ratio Xs (see text for definition) as determined in elementary 
and heavy ion collisions from fits of multiplicities to the statistical hadronization model. 



2.8 Thermalization: How Is It Achieved? 

After discussing the success of the SHM in reproducing particle multiplicities and the 
intriguing universality of its main parameter, the temperature, one is obviously led to 
the question how this can come about. A classical process of thermalization through 
binary collisions between formed hadrons, advocated in heavy ion collisions [35], is 
ruled out in elementary collisions because the expansion rate is fast and hadrons are 
not interacting for a long enough time for this to take place. But even in heavy ion 
collisions peculiar features of the data cannot be explained in a hadronic kinetic picture 
[36] without invoking the predominance of multi-body collisions; since, in this case, 
the hadronic mean free path is comparable or smaller than Compton wavelengths, the 
collisional picture breaks down naturally. 

There is evidence that thermalization occurs at a relatively early stage over a large 
region (i.e. clusters several femtometers wide) in heavy ion collisions, whereas it is a 
late phenomenon (i.e. very close to hadronization) occurring over small (of the order of 
1 fm) distances in elementary collisions. Yet, the agreement between model and data 
is surprisingly accurate in elementary collisions, even more accurate than in heavy ion 
collisions, the only difference being in the level of strangeness phase space saturation. 
Somehow, the hadrons must be born into equilibrium as Hagedorn first pointed out [37] 
and was reaffirmed by others [38j [39] . 

The idea that this thermal-like behavior is of genuine quantum-mechanical origin 
and not related to semi-classical collision processes, is shared by many [32] and is es- 
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Figure 11: Canonical enhancement (defined here as the ratio between the chemical 
factor Z(S — Sj)/Z(S) and its value at some fixed volume Vo) as a function of volume 
for hyperons (from ref. [29]). 



poused in this review. A different point of view was presented in ref. jlQ] where it was 
argued that the thermal behavior could just be mimicked by a matrix element which is 
weakly dependent on the final kinematic variables in ( 1T6]) . a scenario called "phase space 
dominance" . But even this scenario requires stringent conditions on the dependence of 
cluster decay rates on the channel multiplicity (essentially like A N ) [10] such that the 
exponential dependence of production rates on mass is not spoiled. Hence phase-space 
dominance is not less trivial in any way. A possible path to distinguish between the two 
scenarios is provided by the analysis of exclusive rates at low energy, although it must be 
pointed out that the observed identical particle correlations already favors SHM which 
is endowed with a built-in spacial extension, unlike phase-space dominance. 

However, whether it is a proper thermal-statistical equilibrium in a finite volume 
or rather a phase space dominance effect, there must be a profound reason behind this 
phenomenon, which ultimately has to be related to the nature of QCD as a theory with 
strong coupling at large distances. Also, we believe that the intriguing universality of 
the temperature found in elementary collisions as well as heavy ion collisions and its 
resemblance of the QCD critical temperature is not accidental. 

If we assume that post-hadronization collisions are unable to restore equilibrium, 
how can a quantum evolution process ensure it? Several years ago it was pointed out 
that a closed quantum system whose classical counterpart is chaotic and ergodic can 
give rise to thermal distributions provided that the so-called Berry conjecture applies 
|41j . Berry's conjecture [12] essentially states that the high-lying eigenfunction am- 
plitudes ^>(x) in configuration space appear to be random Gaussian numbers and, as 
a consequence, momentum space distribution is microcanonical [13]. This "quantum 
thermalization" mechanism has been invoked to explain the observed thermal-like dis- 
tributions in hadronic processes [H]. Of course, this argument requires that classical 
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QCD is chaotic (as it has indeed been advocated [15]), that Berry's conjecture holds for 
quantum fields, and that it can be applied to a dynamical process like hadronization. 
All of these conditions are non-trivial, but pursuing these ideas further may give rise to 
interesting developments. 

Recently, another appealing idea to explain the universality of thermal features in 
multihadron production has been put forward [3B]- It invokes an analogy between con- 
finement and black hole physics. It is conjectured that the phenomenon of confinement 
is equivalent to the formation of an event horizon for colored signals (quarks and glu- 
ons). Similarly to Hawking- Unruh radiation, the spectrum of hadrons, emitted as the 
result of a high energy collision, is thermal because no information can be conveyed 
from the causally disconnected region beyond event horizon. According to this so-called 
Hawking-Unruh scheme of hadronization, temperature is related to the string tension 
and is thereby universal. Another interesting consequence of this idea is that the extra 
strangeness suppression observed in elementary collisions can be quantitatively explained 
[T7] as an effect of the different strange quark mass. 

These attempts relating the observed thermal features in hadron production pro- 
cesses to quantum chaos or Hawking-Unruh radiation are still in a developmental stage. 
Whether they will keep their promises will be seen in the future. Certainly, both share 
the vision that there is a fundamental quantum mechanical mechanism behind this 
phenomenon and no (or little) room for a classical collisional thermalization process. 

3 Quark Recombination 

The statistical hadronization model provides a successful description of key features of 
hadron production without explicitly invoking the underlying fundamental degrees of 
freedom in QCD, quarks and gluons. Clusters inherit the characteristics of partons from 
which they emerge, but it seems that there is no explicit role for hadronic substructure in 
determining final hadron ratios except for the extra strangeness suppression parameter 
7s which is a direct manifestation of the differences between u, d and s quark masses. On 
the other hand, many hadronization models involve parton degrees of freedom explicitly 
and it seems that some observables in relativistic heavy ion collisions at very high energy 
require the assumption of an underlying parton dynamics. 

In a parton-based approach the goal is to calculate the probability to produce a set 
of hadrons hi, . . ., h n from a given partonic "initial" state C. In general, this problem 
involves interactions of partons at a scale around Aqcd which, as emphasized before, is 
a highly non-perturbative problem and as yet unsolved. In defining this problem, we 
have to overcome another obstacle which is connected to the preparation of the partonic 
"initial" state C. Quarks and gluons are not asymptotic states. Rather, they usually 
appear as intermediate states in the scattering reaction A+B — > C+X — > h±+. . . h n +X'. 
Hence the partonic state C is difficult to "prepare" , and in fact the nature of quantum 
field theory requires us to integrate over all possible states C. In the worst case, C might 
not be well defined at all because the partonic state couples to the other states, A, B, 
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X, etc. in a non-factorizable way. In fact, this will always be true for any state C which 
is not a color singlet, since the final hadron ensemble hi, . . ., h n is colorless and color 
must be exchanged with other parts of the system during hadronization. We call this 
problem the factorization problem. 

Several factorization schemes for hadronization have been developed starting from 
the 1970s, resulting in some successful effective descriptions of hadronization. They 
are usually based on a proof that in a well defined process a simple set of possible 
intermediate parton states C« is giving the leading contribution to hadron production 
in powers of a small expansion parameter. The best known example is single hadron 
production in e + e~or hadronic collisions at a large momentum scale Q. The leading 
contribution to the cross section in terms of powers of 1/Q comes from intermediate 
states with just one single parton C\ = {g}, C2 = {«}, . . . and the cross section o h 
can be factorized into a cross section for producing the intermediate parton Ci and a 
probability D i/h to produce h from C [1711181119] 

<7 fc ~ J]^® D i/h . (27) 

i 

Note that the leading contribution does not involve an interference effect between dif- 
ferent partons in the amplitude and the complex conjugated amplitude and the process 
C{ — > h can be formulated as a probabilistic problem with a probability distribution 
Di/h{z) called a fragmentation function, z is the fraction of the original parton momen- 
tum carried by the hadron h. In general, the functions Di/h(z) can not be calculated 
from first principles (which would be equivalent to fully solving the hadronization prob- 
lem). However, it makes single hadron production at large momentum treatable by 
separating off the hadronization probability for a given hadron which can be measured 
by experiment and is universally applicable to all processes where such a factorization 
holds. 

All hadronization models involving partons have to start from an assumption about 
a well defined parton "initial state". The exact details will depend on the process 
under investigation. Note that the SHM had such a basic axiom as well, defining the 
probability of multihadron states in a cluster. In this section, we discuss the model of 
quark recombination or coalescence. It also involves a plausible but ad hoc assumption 
at the outset. This is less rigorous than the factorization argument that can be given 
for parton fragmentation, but it is phenomenologically very successful. We will focus on 
applications of quark recombination to heavy ion collisions. 

3.1 Parton Fragmentation and Its Limitations 

For a better understanding of recombination it is instructive to investigate situations 
where parton fragmentation is failing. 

Factorization of fragmentation functions off a hard parton scattering cross section 
was first introduced for simple single scale processes like single-inclusive hadron produc- 
tions [SD]- Later the LEP Collider at CERN brought reliable data on hadron production 
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in e + e~ collisions which allowed the extraction of universal sets of fragmentation func- 
tions [5H [52l [53] (see [49] for a modern review). In a physical gauge fragmentation 
functions have a straightforward definition in terms of counting operators a\ l (P)ah{P) 
of asymptotic hadron momentum eigenstates applied to the parton field. E.g., if parton 
i is a quark described by a field ip then we have 

D i/h {z) ~ J dxe- lPx / z (0\^(0)a[(P)a h (P)tP(x)\0) (28) 

where x is a light cone coordinate conjugate to the light cone momentum P. 

There are many situations where factorization involving a vacuum fragmentation 
function could not be proven or is explicitly violated. In some situations it is easy to see 
qualitatively why this is the case. One example is hadron production at very forward 
rapidity in hadronic collisions. A phenomenon known as the leading particle effect [HI] 
can be observed if relative abundances of hadrons at forward rapidity are considered. 
The FNAL E791 collaboration [55] found a large asymmetry between D~ and D + mesons 
in fixed target experiments with 7r~ beams on nuclei. Obviously a fragmentation from 
cor c quarks or from gluons should be nearly flavor blind and the charm quark mass 
sets a (semi-)hard scale that can be used to justify perturbative arguments. However, 
in this case the asymmetry comes from recombination of the c from a cc pair produced 
in the collision with a d valence quark from the beam tt~ remnants. This mechanism is 
enhanced compared to the c+d recombination which involves only a sea quark from the 
7r~ [56]. Obviously the presence of the beam fragments acting as spectators destroy the 
favorable conditions for vacuum fragmentation. 

The leading particle effect also motivated the birth of the coalescence picture. Das 
and Hwa developed a model for recombination of partons streaming forward in hadronic 
collisions [57] . There is no thermalized parton phase in this case, but this simple model 
has many of the features of the recombination models for RHIC we discuss below: 

(1) The input is a (multi-)parton distribution f(pi, . . . ,p n ) whose spectrum and chem- 
ical composition remain unchanged during the hadronization process. 

(2) Partons coalesce into hadrons according to recombination functions $ which play 
the role of squared hadronic wave functions. 

(3) Only valence quarks of the hadrons play a role in the hadronization process. 
The number of recombining hadrons is therefore given by 

JV fc ~$®/(pi,...,20 (29) 

where n = 2 or 3 for mesons and baryons respectively. The first condition is usually 
reinterpreted to mean that the coalescence process is fast enough such that the momen- 
tum distribution and chemical composition of the partons do not react to the depletion 
of partons during recombination. This assumption addresses the question of a well- 
defined "initial" parton state indicated in the last section. The third condition has been 
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the center of many debates. Where are the gluons and sea quarks? While there is no 
calculation of this process from first principles, we do have some qualitative arguments. 
First, the momentum transfer involved is rather small and it is questionable whether 
perturbative partons can be resolved. Second, we know that even in processes that start 
at large momentum scales, degrees of freedom get frozen and give constituent masses to 
quarks [58] . We could therefore argue that a dressing of quarks happens before the final 
recombination step. Third, we have no reason to believe that the hadrons originally 
formed are indeed free hadron states. There will necessarily be a formation time before 
we can actually interpret them as asymptotic states. The last two arguments are related 
to the role of chiral symmetry during hadronization. Chiral symmetry breaking has not 
been explicitly incorporated in recombination models thus far. 

3.2 The Recombination Formalism 
3.2.1 Basic Theory 

Coalescence or recombination of particles appears in a wide array of systems in atomic, 
molecular and plasma physics. As a first approximation the details of the dynamical 
process are usually ignored. Rather, the adiabatic approximation of a projection of the 
initial multiparticle state onto the final coalesced state is considered. This instantaneous 
approximation is widely used in the literature for the case of partons coalescing into 
hadrons. The formalism introduced here is also related to the successful coalescence 
model for nucleons [59j EO] • After first applications to forward hadron production and 
the leading particle effect the coalescence concept was soon applied to heavy ion collisions 

[SDElEaiMl. 

The number of hadrons h coalescing from a partonic system characterized by a 
density matrix p is given by 




Instantaneous here means that the states are defined on a hypersurface which is typically 
either taken to be at constant time, t = const., or on the light-cone t = ±z. Note that 
the notion of instantaneous recombination guarantees that condition (1) from section 
13.11 is automatically fulfilled. 

The instantaneous projection formalism has the conceptual disadvantage that only 
three components of the four-momentum are conserved in the underlying 2 — > 1 or 
3 — > 1 coalescence process. Qualitatively this can be corrected by assuming that the 
participants can scatter off the surrounding particles which can build up or dissipate 
their virtuality. At least in a equilibrium state this should not change the momentum 
distribution or chemical composition of the particles. However, no explicit formalism 
has been developed to include this effect quantitatively. 

More dynamical approaches beyond the instantaneous projection approximation 
have been considered as well [65]. They conserve 4- momentum by allowing a finite 
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width for hadrons and will be discussed later. We first focus on instantaneous projec- 
tion models that were discussed for heavy ion collisions first by Greco, Ko and Levai 
[GKL] [661 EZ1; Fries, Miiller, Nonaka and Bass [FMNB] [681 [70]; Hwa and Yang 
[HY] [7JE1 and Rapp and Shuryak [RS] [71]. 

^From Eq. ( 1291) we can derive an expression for the number of mesons with a certain 
momentum P from recombination [69] 

dN M ^ f d 3 Rd 3 qd 3 r Tjr / r P „ r P \ . . 

Here M denotes the meson and a, b are its coalescing valence partons. W a b and $m 
are the Wigner functions of the partons and the meson respectively, P and R are the 
momentum and spatial coordinate of the meson, and q and r are related to the relative 
momentum and position of the quarks. The sum runs over all possible combinations 
of quantum numbers of valence quarks in the meson, which is usually replaced by a 
degeneracy factor Cm- Eq. (|3T|) is very intuitive for a hadronization hypersurface at 
constant time. It can be made Lorentz-covariant to allow for relativistic kinematics. 

The corresponding formula for baryons, containing 3 valence quarks, is completely 
analogous [69]. Note that Eq. (13T1) implements principle (3) from Section I3TT1 by taking 
into account only the lowest Fock state of the meson. It has been attempted to generalize 
Eq. fl3TT) to include more partons which would be gluons or pairs of sea quarks, accounting 
for the next terms in a Fock expansion [73]. There is no difficulty in including these in 
the model, but data on elliptic flow scaling in heavy ion collisions puts bounds on the 
size of contributions from higher Fock states. 

Let us recall the definition of the Wigner function for a meson consisting of two 
quarks 

$ M (r, q) = J dhe~ is -\ M (r + |) <p* M (r - |) (32) 

where the 2-quark meson wave function in position space ifM can be represented as 

(r i; r 2 |M; P) = e^ p < ri+r ^ M (n - r 2 ) (33) 

The Wigner function of the partons can be defined in a similar way from the density 
matrix p [69] . For further evaluation the Wigner functions have to be specified. Usually, 
the parton Wigner function W a b is replaced by its classical limit, the multi-parton phase 
space distribution f(x\,pi] . . . ; x n , p n ) on the hypersurface of hadronization. The Wigner 
functions of the hadrons are not well known in this particular case. In general the 
wave function of a hadron depends in a non-trivial way on the frame of reference, the 
particular process, and the resolution scale. In heavy ion collisions the resolution scale 
is non-perturbative, and of the order of the temperature at hadronization, T c < 200 
MeV. Therefore the hadron Wigner functions are often modeled ad hoc. Luckily the 
dependence of observables in heavy ion collisions on the shape of the wave function is 
suppressed as we will discuss below. 
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The different groups mentioned above using instantaneous recombination mainly 
differ in details of how they apply the basic equation (13T|) to nucleus- nucleus and hadron- 
nucleus collisions [75]. We will discuss some of their unique features and refer the reader 
to the original literature and reviews [751 ES] for further details. 

3.2.2 Different Implementations of Instantaneous Recombination 

In the implementation by Greco, Ko and Levai [GKL] [66l [67] the full overlap integral 
in Eq. fl3~T|) over both relative position and momentum of the partons is calculated with 
Monte-Carlo techniques. In their case the Wigner function of the hadron is a simple 
product of spheres in position and momentum space with radii correlated through the 
uncertainty relation. 

In the implementations by Fries, Miiller, Nonaka and Bass [FMNB] [68l l69l 170] ; 
Hwa and Yang [HY] [7TJ [72] and Rapp and Shuryak [RS] [7J] Eq. (ED) is simplified by 
integrating out the information about position space. This leads to a pure momentum 
space formalism in which the information about the hadron is further compressed into 
a squared (momentum space) wave function. In a frame of reference in which one 
momentum component of the hadron is much larger than the others the two momentum 
components transverse to the hadron momentum can be integrated out as well, leading 
to a simple 1-dimensional overlap integral. Such a formalism is appropriate, e.g., for 
hadron production at large transverse momentum in the laboratory frame. It is also 
equivalent to the original formulation by Das and Hwa using recombination functions 
[57] . The yield of mesons with momentum P can then be expressed as 



where da is the hypersurface of hadronization. In many cases the emission integral 
over the hypersurface is not calculated explicitly, but replaced by a normalization factor 
proportional to the volume of the hadronization hypersurface. x\ and X2 are the fractions 
of momentum carried by the two valence quarks. 

For the lowest Fock state of a meson the squared wave function is mostly parame- 
terized as [HS] 



The powers on determine the shape while the constant B normalizes the integral over 
$m to unity. Light cone distribution amplitudes suggest = 2 for a light quark 
system [771 ESI [79]. Values a c = 5, a u ^ = 1 are used for the charm-light quark system 
in D mesons in Ref. [H]. The extreme case of the two quarks exactly sharing the 
momentum, $m(xi,X2) = S(x\ — l/2)5(x2 ~~ 1/2), is often considered for schematic 
estimates. Formally this corresponds to oii — > oo with the ratio of the fixed. 

3.2.3 Dynamical Recombination 

One of the drawbacks of the instantaneous approximation is the violation of energy con- 
servation which has to be restored, on average, by interactions in the medium which are 
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not further specified. One early attempt to improve the situation was the introduction 
of a mass distribution for the quarks [80J. The mass distribution can be thought of as 
an effective way to incorporate some in-medium effects. This approach allows to enforce 
both momentum and energy conservation and one finds fairly good agreement with data 
for transverse momentum spectra in heavy ion collisions. 

A promising new direction has recently been taken by Ravagli and Rapp (RR) [651 
IHT] . They advocate a dynamical approach in which the instantaneous projection of quark 
states onto hadron states is replaced by a transport picture which involves ensembles f a 
and fu of quarks and mesons evolving with time. To determine properties of the meson 
distribution starting from given quark distributions one solves the Boltzmann equation 

?-f M (t,P) = -^-f M (t,P) + J d ~j0^ J d 3 xf a (x, Pl )f b (x,p 2 ) 

x a{s)v(p u p 2 )S {3) (P - pi - pa) (36) 

with a gain term describing the scattering of a quark- ant iquark pair a-b into a meson 
state M with finite width T and a corresponding loss term describing the decay of M 
into a pair a-b. Here a(s) is the resonant cross section in Breit-Wigner form with width 
T, v is the relative velocity of the pair and 7p = a/P 2 — m 2 /m is the gamma factor of 
the meson with mass m. 

This implementation naturally conserves 4-momentum and permits, e.g., the study 
of kinetic energy scaling [821 IHT] . The problem of baryon formation is slightly more 
difficult in a transport approach and has not yet been addressed. We should also note 
that even for a finite meson width the creation of light Goldstone bosons, in particular 
pions, is still problematic. They are much lighter than the sum of the participating 
quark masses which are usually taken to be constituent quark masses. 

3.3 Connection with other Hadronization Models 

The recombination model was conceived as a correction to vacuum fragmentation in 
order to include multi-parton effects. Indeed, for a first, intuitive understanding it is 
helpful to view parton fragmentation and parton recombination as two limiting cases 
of hadronization. The former fixes the parton state C to be a single parton in the 
vacuum. The latter assumes the parton state to be a subset of a populated, maybe even 
thermalized phase space, such that valence partons are readily available. This simple 
picture allows a basic understanding of most qualitative features of both mechanisms. 

Recombination is the most effective way to form hadrons from a thermalized parton 
phase space, at least for not too small momentum P. This is easy to see in the instan- 
taneous momentum space formalism, but it can be shown to hold numerically in more 
general cases as well. Let us assume that the partons at moderate to large values of p 
(large enough that the energy E can be approximated by E ~ p) are distributed ac- 
cording to a Boltzmann distribution f a = e~ p ^ T . Fragmentation then leads to a hadron 
spectrum which has an effective slope (z)T with (z) < 1 being the mean value of the 
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momentum fraction z of the hadron from a fragmentation function D(z). Hence the 
resulting slope is steeper than 1/T. On the other hand, hadrons recombining according 
to Eq. (1531) have a distribution 



for mesons, and similar for baryons. Hence we conclude that the recombination formal- 
ism is indeed more effective for large P as long as the spectrum is exponential. 

We also notice that for large momenta (P » M) we recover the results we obtain 
from the (momentum dependent) statistical model. At least this is true if the same set of 
resonances is included. In practice, this has not yet been implemented for recombination. 
Only a few studies have been dealing with resonances in a limited way [83J, E7J IM] . 
Of course, one can argue that the contributions from resonance decay to stable hadron 
spectra are formally vanishing in the limit P — > oo. With this caveat the correspondence 
between recombination and statistical hadronization would be precise. 

Unfortunately, the limit P — > oo is an academic one, since in any realistic system the 
parton spectrum will eventually become a power law f a ~ p~ a for increasing momen- 
tum p. Once this kinematic regime is reached fragmentation becomes more effective, 
producing hadrons with a slightly shifted power law ~ p-a-s^ w } iei[e small expo- 
nent 5 approximates the changed slope after applying the fragmentation function. On 
the other hand the power is doubled or tripled for recombining mesons (~ p~ 2a ) or 
baryons (~ P~ 3a ), respectively. Hence, in any realistic system fragmentation dominates 
at asymptotically large momentum, in accordance with QCD factorization theorems. 

One might ask whether recombination is interesting to consider given its apparently 
small region of applicability. Curiously, in heavy ion collisions at RHIC we find a kine- 
matic region at intermediate transverse momentum between 1.5 and 6 GeV/c where 
recombination still dominates over fragmentation, and the momenta are large enough 
compared to the masses involved such that the simple instantaneous projection formal- 
ism, even without including resonances, gives very reasonable results for hadron spectra 
and elliptic flow. Moreover, unlike statistical hadronization, recombination can make 
certain predictions connecting the kinematic behavior of hadrons and quarks. Ellip- 
tic flow scaling in heavy ion collisions is a famous example. It must be connected to 
hadronization because it exhibits a dependence on the number of valence quarks. We 
will offer a survey of recent experimental results supporting quark recombination further 
below. 

Let us summarize this subsection so far. We have found that there is a formal 
connection between parton recombination and the statistical hadronization model at 
large momenta. Without doubt a successful implementation of parton recombination 
at low momenta in heavy ion collisions should reproduce the statistical model in that 
kinematic region as well. Such a description is not available yet. It would have to 
include energy and momentum conservation, at least in a statistical sense, as well as 
the full range of resonances and their decays. This issue of resonances is also related 
to entropy conservation. Direct coalescence into stable hadrons seems to reduce the 
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number of particles by about a factor two, raising questions about the second law of 
thermodynamics. At intermediate and large momenta, where the current formalism is 
applicable, only a minute fraction of the total number of particles is residing, so that 
entropy is not of any grave concern here. 

We conclude this subsection by discussing some further connections between recom- 
bination and fragmentation. We return to the original question of preparing a well- 
defined partonic state C. The fragmentation formalism puts a cut between long- and 
short- distance physics which is well motivated because we can not calculate the former 
perturbatively, while we can compute the latter. The single parton state is treated like 
an observed final state of the perturbative calculation. However we know that many 
processes during fragmentation are still partonic in nature, not hadronic. Hence we are 
tempted to ask whether we could define a partonic state C later in the development of 
the jet parton shower. A jet shower might include a dozen or more partons which at 
some point feel the pressure from the QCD vacuum to bind into hadrons. This binding 
might well be described again by recombination. 

It is not at all clear that such a parton ensemble in a jet "just before" hadronization 
can always be well defined. However, the cluster picture of statistical hadronization 
seems to support this idea. One can try to choose this as a starting assumption. This 
picture of the fragmentation process has been implemented in the recombination model 
of Hwa and Yang [85]. Instead of fragmenting hard partons directly, they define the 
parton shower of a jet (initiated by a hard parton) through so-called shower distribu- 
tions. These are given by non-perturbative splitting functions Si/j(z) which describe the 
probability to find a parton of flavor j with momentum fraction z in a jet originating 
from a hard parton i. The parton content of a single jet is then allowed to recombine 
and should match the corresponding fragmentation functions. E.g. for a meson M with 
valence quarks a and b 



Di/ M {z) = / dx a dx b 



Si/ a (x a )S i/b [ % ^ ) + (a <-> 6) 



$M{x a /z,x b /z) . (38) 



The shower distributions are not determined from first principles but are fitted to re- 
produce the known fragmentation functions for pions, protons and kaons [85] . 

The power of this approach lies in the fact that the fragmentation part of the hadron 
spectrum is computed with the same recombination formalism that is applied to the 
thermalized part of phase space. It is then very natural to allow coalescence of shower 
partons with thermal or soft partons [72] by defining a total parton distribution 



fa( P ) = fi*(p) + /f owcr (p) =: T(p) + S(p). (39) 

which is subject to the recombination formalism. For mesons we obtain two known 
contributions. The term ~ SS corresponds to the coalescence of two partons from the 
same jet and reproduces fragmentation, the term ~ TT recombines two soft partons and 
reproduces the soft recombination spectrum discussed in the previous subsection. The 
mixed term ~ TS is new and describes the pick-up of a soft parton by a jet parton to 
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form a meson. For baryons two of these mixed terms, ~ TTS and ~ TSS, exist. These 
soft-hard coalescence terms make important contributions in some cases. They can, e.g., 
explain the final state dependence of the Cronin effect observed in d + A collisions at 
RHIC [861 [87], and lead to jet-like correlations between hadrons [88| 1891 190| I9T] . 

To summarize, there is reason to believe that parton coalescence could be an over- 
arching qualitative concept for hadronization, based on its phenomenological success. 
However, so far it is only defined in a naive parton model sense. Beyond that it is 
not clear that the parton content of an arbitrary scattering of hadrons or nuclei before 
hadronization can be written down in a well-defined way. E.g. in Eq. fl39|) . it is not 
clear why hadrons forming at very different times in the soft and hard sector can be 
treated together. It is also not completely understood how the scale dependence of the 
fragmentation functions can be interpreted correctly [92J. 

3.4 Experimental Evidence for Recombination 

We have already discussed the leading particle effect and mentioned the Cronin effect in 
hadron-nucleus collisions. Here, we want to review some of the support for recombination 
from heavy ion collisions. Three early results from the Relativistic Heavy Ion Collider 
led to the conclusion that hadrons at intermediate transverse momentum px (1.5 GeV/c 
< pt < 4 . . . 6 GeV/c), are produced by recombination of quarks. 

(a) The baryon-to-meson ratios in central Au-Au collisions was found to be enhanced [931 
[9l]. The proton/pion ratio ~ 1 is incompatible with expectations from fragmen- 
tation. 

(b) The nuclear modification factors Raa and Rcp, i-e. the ratio of yields in central 
Au-Au collisions compared to pp and peripheral collisions, resp., scaled by the 
number of binary nucleon-nucleon collisions, is not suppressed for baryons [931195]. 

(c) The anisotropy of particle production in azimuth relative to the reaction plane, 
called the elliptic flow parameter f 2) scales in a universal way for mesons and 
baryons [HS1 Ml E3 EE] ■ 

3.4.1 Hadron Spectra and Ratios 

Fig. [T2l shows the measured anti-proton/pion [93] and A/ Kg [94] ratios as functions of 
transverse momentum for various centralities and collision systems. At intermediate 
Pt the baryon-to-meson ratios in central Au-Au collisions are significantly larger than 
those in e + e~ [99] or pp collisions [100] . These results were the first clear indication 
that hadronization proceeds differently in Au-Au collisions and pp collisions even in a 
regime traditionally expected to be dominated by jet fragmentation. 

Fig. [13] shows anti-proton/pion (left) and A/K s (right) ratios compared with instan- 
taneous coalescence calculations from Greco, Ko and Levai [GKL] [671 HOlj and Fries, 
Miiller, Nonaka and Bass [FMNB] [691 1102] . Despite some deviations these calculations 
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Figure 12: Left: p/vr~ ratios measured in central Au-Au collisions at ^s NN = 200 GeV 
at RHIC, compared to measurements from e + e~and pp collisions. Right: The ratio 
A/2K$ for central and mid-central Au-Au collisions at y/s NN = 200 GeV measured by 
STAR. The pj ir~ ratio from pp collisions from STAR is shown for comparison. 

reproduce the characteristic peak structure very well. This peak comes from the inter- 
play of soft production which leads to steeply rising ratios at low p?, and low ratios at 
large p? coming from jet fragmentation. Recombination gives an explanation why soft 
production dominates out to intermediate px- 

The different behavior of baryons and mesons is also expressed in different nuclear 
suppression factors. Fig. [TD shows the nuclear modification factor Rcp measured at 
RHIC for various identified hadron species. If a nuclear collision is a simple superposition 
of binary nucleon collisions, Rcp will equal one. Therefore the suppression of Rcp at 
large pt is a signature for jet quenching pL03j . The baryons (A + A, S + H, and f2 + Q) 
[941 1104j . show significantly less suppression than mesons (kaons or 0) [105[ 1106] ). The 
same behavior was found for protons and pions |107j . The second important observation 
is that the behavior is universal for all baryons and all meson separately, and largely 
independent of mass. Heavy <fi mesons behave like pions, not like the equally heavy 
protons or A baryons [1061 [82] . This is a clear indication that the decisive factor here is 
the number of valence quarks. 

Fig. [15] shows spectra from the FMNB coalescence model [69] for neutral pions, 
kaons, protons and hyperons in central Au-Au collisions compared to data from RHIC 
[10411107] . The recombination region is visible as an exponential slope at intermediate pr 
with a transition to a power-law shape at higher pr- Note that the FMNB model simply 
adds recombination and fragmentation contributions. Nevertheless the description of the 
data is quite good above 2 GeV/c. Note that the transition from an exponential shape 
to a power-law shape happens at a higher pr for baryons than it does for mesons. This 
is expected and due to the inherent suppression of baryons in vacuum fragmentation 
which makes coalescence competitive up to quite large px- 
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Figure 13: Ratios of baryon yields to meson yields for central Au-Au collisions at 200 
GeV. The GKL and FMNB calculations for p/tt~ (left) and A/2K® are compared to 
STAR and PHENIX data. 



3.4.2 Elliptic Flow 

Recombination makes a straightforward prediction for the elliptic flow of hadrons 
[108] . using the elliptic flow of quarks as input. This behavior can be easily derived in 
the instantaneous momentum-space formalism. Let us assume that the elliptic flow of 
a set of partons a just before hadronization is given by an anisotropy v 2 (pt) at mid- 
rapidity. The phase space distribution of partons a can then be written in terms of the 
azimuthal angle as 

/a(pr) = Upt) (1 + 2v a 2 (p T ) cos 20) , (40) 

where odd harmonics are vanishing due to the symmetry of the system and higher 
harmonics are neglected. / is the distribution averaged over the azimuthal angle 0. 

For a meson with two valence partons a and b and for small elliptic flow ^ < 1 one 
has [69] 

Mi \ J d(f)cos(2(f))dN M /d 2 p T 
J dcpdN M I d z p T 

dx a dx b $ M (xa,Xb) [v%{x a pr) + v 2 (x b p T )] ■ 

In the case of a very narrow wave function in momentum space (a — > oo) this leads to 
the expression 

<M = v a 2 (x a p T ) + v b 2 {x b p T ) . (42) 

with fixed momentum fractions x a and x b (x a + x b — 1). 

Thus for hadrons consisting of light quarks which exhibit the same elliptic flow before 
hadronization we arrive at a simple scaling law with the number of valence quarks n: 

v h(p T ) = nv a 2 (p T /n) . (43) 
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Figure 14: Nuclear modification factors (Rcp) for various identified particles measured 
in Au-Au collisions at y/s NN = 200 GeV by the STAR collaboration. The values for 
baryons and mesons fall into two separate bands (indicated by lines to guide the eye) 
with the baryon Rcp larger than the meson Rcp- 



This scaling law had first been found in data from RHIC and was quickly interpreted 
as a coalescence artefact |109j . Since then this connection has been solidified in al- 
most all coalescence models [691 H10[ 1111] • Eq. (}4"2"|) has also been used to extract the 
elliptic flow of heavy quarks from measurements of D mesons |112l 1113] . The line of 
thought presented for elliptic flow can be extended to harmonics beyond the second 
order. Generalized scaling laws for the 4th and 6th order harmonics have been derived 
in Ref. [TT4] . 

Fig. [16] shows data on v 2 scaled by the number n of valence quarks in a given 
hadron as a function of pr/n for several species of identified hadrons at y/s NN = 200 
GeV |106[lll5j . To investigate the quality of agreement, the data from the top panel are 
scaled by a polynomial fitted to the universal curve and plotted in the bottom panel. At 
low pt hydrodynamics |116l 11171 1118] predicts an ordering of v 2 by hadron mass which 
is confirmed by data [951 E3 1119] • That is the reason why valence quark scaling of pt 
is not working for pr/n < 0.5 GeV/c. The two phenomena can be reconciled by first 
scaling the vi results in transverse kinetic energy Et = mr — mo, and then applying a 
scaling in n [H2] • The result is also shown in Fig. []Jjl The scaling at low — m holds 
to high accuracy. 

The experimentally observed scaling is a spectacular success for the recombination 
picture. However, it is also the source for some open questions. First, the derivation 
of the scaling law in the instantaneous momentum space formalism is rather simplistic 
and neglects both energy conservation and space-momentum correlations in the quark 
distributions. The observed scaling seems to put tight constraints on the latter |120[ 
12 lj. The instantaneous formalism is also ill-prepared to distinguish between valence 
quark scaling in px and Et because of the partial violation of energy and momentum 
conservation. Recent work using dynamical recombination with energy conservation 
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Figure 15: Hadron p^-spectra at midrapidity from 200 GeV central Au-Au collisions. 
The curves show the recombination and fragmentation components of the spectra ob- 
tained in the FMNB formalism along with the total which compares well with the data. 

and realistic space-momentum correlations found that recombination respects scaling 
in Et [HI]- Future work in this direction is required to understand all aspects of the 
scaling phenomenon and its breaking |110[ I112[ [73| 1122] . It will also be interesting to 
see how prominent recombination effects will be at the higher energies reached at the 
Large Hadron Collider LHC [T23] . 

4 Summary and conclusions 

In this chapter, we have presented two models of hadronization motivated by phe- 
nomenological observations in elementary and nuclear collisions. 

The statistical model postulates that hadronization proceeds through the formation 
of massive colorless objects called clusters or fireballs, according to the general ideas of 
preconfinement of other cluster hadronization models [5]. They are formed in fragmen- 
tation in e + e~ collisions just as they are in relativistic heavy ion collisions. While the size 
and charge distribution of these clusters differ between collision systems, each cluster 
gives rise to multihadronic states in a purely statistical fashion. Mathematically this can 
be described by a density matrix of localized hadronic states projected onto eigenstates 
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Figure 16: Top left panel: The elliptic anisotropy parameter v?, scaled by quark number 
n and plotted vs pr/n. A polynomial curve is fit to all the data. The ratio of Vijn 
to the fit function is shown in the bottom left panel. Right panels: same but with pt 
replaced by ijit — m 

with all possible hadronic quantum numbers. Interactions are, to a first approximation 
accommodated by including hadronic resonances (ideal hadron-resonance gas model). 

The statistical model has been applied to a wide variety of small and large systems 
which differ by more than two orders of magnitude in the collision center of mass energy 
y/s, revealing intriguing universal features of the hadronization process. The assump- 
tions of the model were found to hold in a remarkable way for relative abundances and 
transverse momentum spectra of both light and heavy flavored species. For large enough 
ensembles of clusters a canonical statistical description can be used, which allows the 
introduction of a temperature T of the clusters which is typically between 160 and 170 
MeV, corresponding to an average energy density of e ~ 0.5 GeV/fm 3 in a single cluster. 
An additional strangeness suppression parameter 7 S is needed to reproduce the data, 
which is found to be between 0.5 and 0.7 in elementary collisions, while 7*, ~ 1 in central 
heavy ion collisions at RHIC. The origin of such an additional parameter and its relation 
to the strange quark mass are not clarified yet. 

Clusters of hadrons at hadronization capture several important features of hadron 
production. It remains unspecified in this model how to connect it to the previous 
partonic phase in a microscopic picture. This step is complicated by the nature of the 
final partonic stage as an intermediate state in a large quantum system. This problem 
has been successfully resolved only for two simple cases: for single partons in high energy 
collisions described by vacuum fragmentation, and for the thermalized parton phase in 
heavy ion collisions which hadronizes through parton coalescence. 

We have discussed the quark recombination model for heavy ion collisions and we 
have shown that this model recovers important features of the statistical model if applied 
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to a thermalized parton phase. Experimental evidence for quark recombination comes 
from the Relativistic Heavy Ion Collider, in particular the spectacular valence quark 
scaling of elliptic flow. This is an important step to prove the existence of collective 
effects on the level of partons in heavy ion collisions |124j . While instantaneous recom- 
bination cannot capture the bulk of the hadronization process and is only applicable 
for momenta larger than 1 or 2 GeV/c, more recent models based on transport descrip- 
tions are able to address the urgent questions of energy conservation and entropy during 
hadronization. We have also shown a possible way to reconcile parton fragmentation 
with statistical hadronization via a recombination model for the full shower of partons 
created inside a jet cone. 

Quark recombination and the statistical hadronization model offer a successful phe- 
nomenological description of many aspects of hadron production in elementary and 
nuclear collisions. They are based on few core principles which seem to capture key 
universal features in the process of hadron formation. While they do not answer the 
central questions of confinement and chiral symmetry breaking, which are of fundamen- 
tal importance for QCD, they allow us to understand some of the key properties of QCD 
in the non-perturbative regime. 
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